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Abstract. 

In this letter we study chaotic dynamical properties of an asymmetrically coupled 
one-dimensional chain of maps. We discuss the existence of coherent regions in terms 
of the presence of defects along the chain. We find out that temporal chaos is 
instantaneously localized around one single defect and that the tangent vector jumps 
from one defect to another in an apparently random way. We quantitatively measure 
the localization properties by defining an entropy-like function in the space of tangent 
vectors. 
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Spatiotemporal chaos has been recently studied in many fields, as Benard 
convection, optical turbulence, chemical reaction-diffusion systems, and so on [1] . One 
of the paradigms of spatial chaos is the chaotic evolution of a spatial pattern. Order 
and chaos are often considered to be opposite notions in nature. One may think of 
order as a stable and regular behaviour, and of chaos as unstable and erratic behaviour. 
Then one may be tempted to conclude that spatial order cannot exist in temporally 
chaotic systems. There are, however, many examples where spatial order coexists with 
temporal chaos, e.g. convection, boundary layer, shear flows. 

Generally, it is difficult to study this spatio-temporal behaviour from the full 
equations, e.g. the Navier Stokes equations in fluids. As a consequence, to gain 
more insight and to facilitate computations, model systems with a set of coupled 
dynamical systems, coupled-map lattice, have been introduced. These are a crude, 
but nontrivial, approximation of extended systems with discrete space and time, but 
continuous states [1] . The simplest example is given by a set of N continuous variables 
Xi which evolve in (discrete) time as 

Xi(n + 1) = (1 - ai - pi) f{xi(n)) 

(1) 

+ aifixi-^n)) + pif(xi + i(n)). 

Usually, the coupling constants are assumed to be independent of the site i and equal, 
i.e. ai = Pi = 7, and periodic boundary conditions are taken. 

In some physical problems, e.g., shear flow, boundary layers or convection, there is 
a privileged direction. This can be introduced in the model (1) by taking asymmetric 
couplings [2,3]: 

ot-i = 7i Pi = 72- (2) 

The system (l)-(2) with periodic boundary conditions shows convective instability 
[4]. In a reference frame moving with a constant velocity v in some band [fmin, ^max] 
an initial perturbation 5xk(0) grows esponentially with a rate given by a co-moving 
Lyapunov exponent X(v): 

Sx k+vn (n)^Sx k (0)e x ^ n . (3) 

The situation becomes much more intriguing and interesting if one introduces non- 
periodic boundary conditions. See, for example, Ref. [2] in which an unidirected lattice 
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map is used to mimic convective turbulence. Recently Aranson et al [5] studied the 
system (1) with the following coupling constants: 

a* =71, Pi =72 for % = 2, . . . ,N - 1 

ai =0, /?i =72, qiat = 7i) = 0, (4) 
7i >72 

and /(#) = R — x 2 with = 1.67, for which the behaviour of a single map is chaotic. 
The system (1), (4) has the stable uniform solution, 

Xi(n) = x(n), x(n + 1) = f(x(n)). (5) 

Aranson et al [5] found that starting from a randomly nonuniform initial conditions, 
after some iterations the Xi(n) become partially synchronized: Xi(n) ~ x(n) for % < l c , 
while for i > l c the Xi are spatially irregular. The finite coherence length l c is due 
to the numerical noise. Indeed, they found that l c increases logarithmically with the 
noise level in the numerical calculations. 

This situation seems rather pathological, since the main features of a chaotic 
system should not be destroyed by the presence of a small noise. 

In order to gain intuition on this behaviour, we have studied the evolution of the 
tangent vector of the system (1), (4). The evolution law of the tangent vector z is 
obtained linearizing equation (1) about the trajectory: 

Zi(n + 1) = (1 - ^ - Pi) g{xi{n)) Zi{n) 

(6) 

+ aig{Xi-x[n)) z*_i(n) + Pig(x i+1 (n)) z i+1 (n), 

where g(y) = df(y)/dy, and the coupling constants are given by (4). 

A numerical analysis reveals that if the tangent vector is initially localized in the 
irregular zone of the system, i.e. Zj(0) = for % < l c , then one has a behaviour similar 
to the convective instability. The vector translates, while growing, in the direction of 
increasing i. However, unlike the case of periodic boundary conditions, as soon as the 
perturbation reaches the boundary i = N, its amplitude decreases rapidly to zero. 

On the other hand, if ^(0) 7^ for some i < l c , the instability does not travel and 
there is an exponential growth of z(n) in the synchronized part of the chain. This 
means that the chaotic part of the system is the spatially coherent one. 

A measure of the degree of chaos is given by the Lyapunov exponents. They 
measure the growth of the tangent vector for large time, and hence they are ruled 
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by the synchronized part of the system. As a consequence, since l c depends on the 
noise level in the numerical calculations, we expect that also the Lyapunov exponents 
should be affected. This is confirmed by the numerical analysis. For example, in the 
case of N = 70 we find only one positive Lyapunov exponent, and the others negative, 
in a four-byte precision calculation, while an eight-byte calculation leads to several 
positive exponents. 

As stated above, a similar scenario, is rather pathological in the framework of 
chaotic systems, since one expects a sort of structural stability for the main features. 
In our opinion the origin of this odd behaviour is in the open boundary conditions. 
Indeed, for a one-dimensional system open boundary condition represents a very strong 
"defect" . The chain is broken somewhere. For this reason we have studied a slightly 
modified version of the model which in some sense interpolate between the periodic 
and the open boundary conditions. 

We consider the system (l)-(2), with 71 > 72 and periodic boundary conditions, 
and we introduce some "defects" by changing the coupling constants at certain points 
i = ki, ■ ■ ■ , /cm (with M <C N) of the lattice. 

We have considered the following class of defects: 

a; = 7i 5 Pi = 72, but «i=72, Pi = 7i) if i = kj. (7) 

The dynamic evolution of a defects, i.e. of a variable with the exchanged couplings, 
Xfc . , is essentially given by the single map solution x, equation (5), as can be seen in a 
return map, see Fig. 1. Thus the motion of the defects are not very much influenced 
by the other variables. Yet there is a small synchronization length, / c , of the variables 
near the defects: Xi(n) ~ Xfc. (n) for i — kj < l c . Unlike the case discussed above, 
this synchronization length does not depend on the noise in the numerical calculation. 
The doubling of the precision in the computation does not produce any change on l c . 

In spite of this apparently simple behaviour, the dynamics is very interesting if one 
studies the evolution in the tangent space. After a short transient, independent on 
both £c(0) and z(0), the tangent vector z(n) get localized around one of the defects. 
It stays there for a while, and then 'jumps' on one of the two near neighbour defects. 
It stays there for a while, and then jumps again, and so on. We have, therefore, a 
sort of random walk of the position of the tangent vector among the defects. In the 
meantime its modulus grows exponentially revealing that the motion is chaotic. 
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One can generalize the analysis by considering the full set of Lyapunov vectors 
z^\n), with j = 1,...,N, associated with the Lyapunov exponents Xj. Each 
vector z^\n) evolves according to (6), but with orthogonal initial conditions, i.e. 
z^(0) _L z^ \0) if j 7^ j'. The Lyapunov exponents measure the exponential growth 
of the volume individuated by the Lyapunov vectors. More explicitly, 

Al + A2 + " ' + K ~ n ln |zW(0) A z^HO) A • ■ • A z(*>(0)| ' (8) 

where A indicates the external product. 

The study of the time evolution of the tangent vectors for the system (1), (7) 
reveals that z^\n), for j < M has the same qualitative behaviour of z^'{n), while 
for j > M does not localize. Figure 2 shows (|2 r l - J ' ) | 2 ) as a function of i, where 
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A quantitative characterization of the degree of localization may be obtained from 
the entropy [6] 

N 

^) = _^/|^)|2\ ln /|^)|2). (10) 



i=l 

We find that when j < M, depends on the density of the defects, p = M/N, 

and that there is a critical density p c ~l c /N. For p < p c we have: 

H j -\nN~ clnp, (11) 

where c depends on 71 and 72, while for p > p c all the entropies IP saturate to a 
constant value, see Fig. 3. 

Finally we note that the Lyapunov exponents tend to cluster according to the 
number of the defects, so that the spectrum Aj depends only on the density of the 
defects, and exhibits the following behaviour 

Xi ~ Xxhp (i/N) , (12) 

where Ai depends only on the parameter of the local map and h p (i/N) is a p-dependent 
function of i/N, see Fig. 4. Formula (12) is well known [7] in the case of coupled maps 
without asymmetric coupling and in the limit of large N. 
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Also these results do not depend on the noise level in the numerical calculations. 

From our analysis, by means of Lyapunov exponents and tangent vectors, we have 
the following scenario. The system is driven by the defects, which practically are not 
influenced by the rest of the system; moreover chaos is concentrated on these defects, 
and the density of the defects is the relevant property to determine the features of 
the Lyapunov exponents. Something similar to a deconfining transition concerning 
the localization properties of chaos appears only for density of defects large enough to 
permit considerable overlaps between the coherent regions of two adjacent defects. 

We conclude by noting that the main features of this scenario hold also for other 
classes of defects, e.g.: 

«i = 7i) Pi = 72, but a* = Pi = if i = kj. 

This means that as long as the chain is not broken (like in the case of Aranson et al.) 
our results are quite robusts and everything is independent of numerical precision. 
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Figure captions 

Figure 1. Xi(n + 1) vs Xi(n) for i = 48, 50, 52, 70 when N = 100 and the defect is 
localized at i = 50. 71 = 0.7, 72 = 0.01. 

Figure 2. (\z^ j) \ 2 ^ vs i for j = 1 (a), j = 2 (6) and j = 3 (c); AT = 300 and two 
defects are present at i = 70 and i = 220. The evolution time is T = 50000. 

Figure 3. - In iV vs p for iV = 200, 300, 500, 800, 1000 

Figure 4. Higher positive Lyapunov exponents, A;, vs i/N in the cases: iV = 100 
(cross), iV = 150 (diamond) and iV = 200 (square), with a constant density of defects 
p= 1/50. 
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